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We present the first numerical QCD bound-state calculation based on a renormal- 
ization group-improved light-front Hamiltonian formalism. The QCD Hamiltonian 
is determined to second order in the coupling, and it includes two-body confin- 
ing interactions. We make a momentum expansion, obtaining an equal- time-like 
Schrodinger equation. This is solved for quark-antiquark constituent states, and 
we obtain a set of self-consistent parameters by fitting B mesons. Applying the 
approach to charmonium spectra leads to a prediction of the hyperfine splitting 
between J/ip(lS) and r/ c (lS) which is in good agreement with experiment. 



1 Introduction 

Recently, a new approach to QCD bound states has been proposed |], ||. The goal is 
to build a bridge between QCD and a constituent quark model (CQM). It has been 
argued that it is convenient to use a light-front formulation of the theory, because on 
the light-front it is possible to make the vacuum trivial simply by implementing a small 
longitudinal cutoff. As a result, all partons in a hadronic state are connected to the 
hadron, instead of being disconnected excitations in a complicated medium. The price 
to pay is a considerably more complicated renormalization problem. 

The new approach consists of two (major) steps. The first step is to find an effective 
Hamiltonian, starting from the canonical light-front Hamiltonian regulated by a large 
cutoff (for details on the cutoff scheme we refer the reader to refs. ||). Unitary 
transformations are used to find counterterms which remove the dependence on the large 
cutoff, and to bring the Hamiltonian towards band-diagonal form with respect to free 
light-front energies. These transformations form a renormalization (semi)group. This 
is repeated until the Hamiltonian is band-diagonal with respect to a typical hadronic 
scale. At the end of this first step the Hamiltonian is still a complicated field theory 
Hamiltonian, but it does not couple states which differ in their free light-front energies 
by more than the hadronic scale. Instead, it contains effective potentials. It was shown 
that the effective potentials contain a Coulomb and a logarithmic confining potential 
already at order g 2 . 

In the second step we want to solve this Hamiltonian and find its spectrum. We 
divide the effective Hamiltonian into a part Hq which is solved nonperturbatively, and 
a part V which is then calculated in bound-state perturbation theory. We want to 
choose Ho so that it is manageable (i.e. something that we can solve) and we want it 
to contain the essential physics. Taking hints from the constituent quark model, we 
include two body potentials and use constituent masses in the Hq, but do not include 
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emission and absorption. Any approximations can always be done by adding a term 
to Hq and subtracting it from V which is treated in bound state perturbation theory. 
Not including the emission and absorption in the Hq has an important consequence: 
different Fock states decouple. We can thus solve few body problems. The consistency 
of this procedure has to be checked in bound state perturbation theory. 

In this talk I present one of the simplest QCD bound state calculation based on 
similarity transformations || and using coupling coherence [||. We find the effective 
Hamiltonian to order g 2 , and then solve for qq bound states. The calculation is carried 
through for quarks of arbitrary but nonzero masses. At the end, we concentrate on 
mesons containing at least one heavy quark. 



2 The effective Hamiltonian to order g 2 

The effective Hamiltonian, which generated by the similarity transformation to order 
g 2 , is band-diagonal in light- front energy with respect to a hadronic scale and it 
can be written as: 

ffeff = tffree + V\ + V 2 + V 2 eff , (1) 

where -fff rec is the light-front kinetic energy (we remind the reader that the light-front 
kinetic energy of a particle with transverse momentum p and longitudinal momentum 

J_ 2 2 

p+ j s p +ra ^ Yy is O(g) emission and absorption with nonzero matrix elements only 

between states with energy difference smaller than the hadronic scale V 2 is 0(g 2 ) 
instantaneous interaction, and V 2 eff includes the effective interactions generated by 
similarity, also 0(g 2 ). The effective interactions generated to this order contain one- 
body and two-body operators. In particular, the effective one-body operator is: 
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where x a = -pq- is the longitudinal fraction of the momentum carried by the constituent 
under consideration, m a is its mass, P + is the total longitudinal momentum of the state, 
V + is the longitudinal scale required in the cutoff by dimensional arguments, and e is 
an infrared cutoff which is to be taken to zero. The divergence in the effective one-body 
operator exactly cancels against the divergence in the effective two-body operator if the 
state is a color singlet. 

The effective two-body operators have the following matrix elements between states 
containing a quark of momentum pi and an antiquark of momentum fcj, i = 1, 2 referring 



2 



to the initial and final state, respectively: 



-g\u{p 2 , er 2 )7 M u(pi, a 1 )v(k 2 , X 2 h u v(k ll Ai)(T Q T b ) 
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where cr, and Ai are light-front helicitics of the quark and antiquark, respectively, u(p, a) 

±2 

and v(k, A) are their spinors, D fllJ (q) = 3 -p2-r] t _ l 7] lJ + -^ 



propagator in light-front gauge, ry M = (0,?7-| 



+ VuQ-^fj) -g^v is the s luon 

1,0,0), q = pi — p2 is the exchanged 



momentum and q 

Do = ko - fcr - q 



D\, D 2 are energy denominators: D\ = p 1 — p 2 — q and 



3 Bound state calculation 

For the purpose of the bound-state calculation we divide the effective Hamiltonian into 
two parts: H which is solved nonperturbatively, and V = H e ff — H which is solved in 
bound-state perturbation theory. In Hq, we include nonrelativistic limit of the kinetic 
energy, the self-energies, and the rotationally symmetric part of the rj^rj^ term of the 
two-body interactions (both instantaneous and generated by similarity) |^]. Let the 
masses of the constituents be m a and m b , and 



M ab = m a + m b . 



(4) 



In the nonrelativistic limit, the light- front scale is naturally replaced by C = ^prjj- 
which carries dimension of mass ||. 
The Hamiltonian Hq is: 
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Hq — 2M a b 
where m is the reduced mass and 
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where 7 is Euler constant. £ contains the finite shift produced by the self-energies after 
subtracting terms needed to make the confining potential vanish at the origin: 
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3.1 Schrodinger equation. 

We now want to find the mass of a qq bound state and its wave function ip(K ± , x): 

cPk 1 - dx 



|p > = / 
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2(2n) 3 y/x(l - x) 
We use a Lorentz-invariant normalization for the states: 

(P'\P) = 2(2ir) 3 P + 6 3 (P - P'), 
and the wave function is normalized to one: 

<,2 ' ;± dx m^,x)\ 2 = i. 



(8) 
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The bound state satisfies: 

H Q \P) =M 2 \P) , (9) 
where M 2 is the invariant mass of the bound state. Let the mass of the bound state be 
M 2 = (m a + m b ) 2 + 2(m a + m b )E , (10) 

which defines E. 

After substituting for H and M 2 in equation (9), some straightforward algebra 
leads to a bound-state equation for the wave function ip: 
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It is convenient to use a dimensionless separation 1Z = Cr that naturally arises in the 
confining piece of the potential, and to absorb — E into a definition of the eigenvalue 
E of the Schrodinger equation. When extracting the bound state mass, — E has to be 
subtracted. The bound-state equation in dimensionless form is: 



+ CaC F (-V conl (U) + V coul (U) 



(12) 



Multiplying both sides of the equation by 2m / C 2 and introducing a dimensionless cou- 
pling and eigenvalue: 



c = 



2maC'F 
C ' 
2mE 



C 2 



we obtain a Schrodinger equation, which depends only on dimensionless variables: 
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ip(Tl) = eip{Tl) . 



(13) 
(14) 

(15) 
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This form is advantageous for numerical study, but moreover, it is quite general - 
one obtains an equation of this form for any quark-antiquark systems and any choice of 
the confining potential in the nonrelativistic limit, regardless of the masses, providing 
they are nonzero. For different systems C, c, e would differ, but the resulting dimen- 
sionless Schrodinger equation will be the same. Thus in the leading order, qualitative 
characteristics of spectra depend only on one particular combination of the masses and 
the coupling, as seen from equations (13) and (14) ||. 

3.2 Results and conclusion 

We refer the reader to ref. (^] for more details on what one can learn from this dimen- 
sionless Schrodinger equation, and the dimensionless results. We applied the approach 
to B mesons 0, and we obtained a set of self-consistent parameters. Applying the 
approach to charmonium, we used IS, IP and 2S levels to fit our parameters (m c = 1.5 
GeV, A = 1.7 GeV, a = 0.5) and then predicted the hyperfine splitting in the charmo- 
nium ground state (the splitting between J/^(1S) and n c (lS) ) to be about 0.13 GeV 

Our study shows that the new approach proposed by Wilson et al. is promising. 
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